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Abstract 

If neutrinos are massive, ‘lepton flavour eigenstates’ are absent from the amplitudes of all 
Standard Model processes. Measurements of r(7T —> en)/r(7r —► fiv) and of the MNS matrix 
are shown to exclude the production of such states in pion decay. Giunti’s objection to this 
conclusion is shown to be invalid. It is pointed out that ‘coherence’ is a property not of produc¬ 
tion processes but of the path amplitudes of quantum experiments. Incoherent production of 
neutrinos of different mass modifies the usual oscillation phase derived assuming production of 
a lepton flavour eigenstate. The new prediction can be tested in long baseline experiments by 
comparison of ‘zq disappearance’ of neutrinos from pion and kaon decay, or by searching for a 
neutrino energy dependence of the phase of atmospheric neutrino oscillations using future data 
from the IceCube neutrino detector. 
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In the Standard Electroweak Model, the coupling of a charged lepton, of generation 
i and a massive neutrino, zq, of generation j to the W-boson is described by the lcptonic 
charged current: 

MCC)“<“ = V * ( 7„(1 - 7(1) 

where U i3 is the Maki-Nakagawa-Sakata (MNS) [1] charged lepton flavour/ neutrino mass 
mixing matrix. Table 1 shows the elements of this matrix as estimated [2] from exper¬ 
imental measurements of atmospheric and solar neutrino oscillations. The non-diagonal 
nature of this matrix gives evidence for strong violation of generation number (or lepton 
flavour) by J ll {CC) levt . Conservation of generation number corresponds to a diagonal 
MNS matrix with zq = zq, z/ 2 = zq and z/ 3 = zq. This is the conventional massless neu¬ 
trino scenario. With massive neutrinos and a non-diagonal MNS matrix, ‘zq’, ‘zq’ and 
‘zq’ do not exist as physical states. That is, they do not appear in the amplitude for 
any physical process. This is because the neutrino charged lepton couplings to the W are 
completely specified by J^[CC) lept which contains only the wavefunctions of the neutrino 
mass eigenstates zq, z/ 2 and z/ 3 . 
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Indeed, the introduction of such ‘lepton flavour eigenstates’ as a linear superposition 
of neutrino mass eigenstates leads to predictions that are excluded by experiment. To 
show this, the pion decays: n —> ev and n —> (iv may be considered. The invariant 
amplitude to produce the anti-neutrino mass eigenstate, V *, in association with a charged 
lepton l~ ( e~ or ji~) is: 


G — 

Mu = -j=f 7T rriT r V ud 'ip e (l - j 5 )U ti ‘ 0 ^ 


( 2 ) 


where G is the Fermi constant, f n and m n are the pion decay constant and mass re¬ 
spectively, and V ud is an element of the CKM [3] quark flavour/mass mixing matrix. 
Introducing a ‘lepton flavour eigenstate’, 0^, as a coherent superposition of the mass 
eigenstates ip Vi : 

= Ua'tpv ! + Ut2i/> V 2 + (3) 

gives, for the invariant amplitude of the decay process 7r“ —> £~Vu 


M( — M^Un + M^ 2 U( 2 + Mf 2 Ui3 


(4) 


where (D is for ‘diagonal’) is given by Eqn(2) with the replacement Uu = 1. If all 
neutrino masses are much smaller than the electron mass, it follows that: 


Mn ~ MYo ~ MY, ~ M 
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where M^ 0 is the invariant amplitude calculated on the assumption of massless neutrinos. 
With two flavour (e-/x) mixing so that: 

Uel U e2 

u„ 1 u ll2 


cos 0 i 2 sin $i 2 
— sin $12 cos $12 


( 6 ) 


and f/ e3 = 3 = 0, Eqns(2)-(6) lead to the prediction: 
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Allowing for radiative corrections [4, 5] the world average experimental value R e /^ = 
(1.230 ± 0.004) x 10~ 4 [6] leads to a determination of the mixing angle $12 via the 


relation: 


The result is: 


or 


1 + tan $12 


= 0.9976 ±0.0032 


1 — tan $12 
tan 2 $12 = (3.6± 3 6 2°) x 10“' 
tan 2 $12 < 3.0 x 10“ 6 at 95% CL 


( 8 ) 


These results are evidently incompatible with the two-flavour mixing results obtained 
from the study of solar neutrino oscillations [2] where $q ~ $i 2 : 

0.22 < tan 2 $Q < 0.71 (LMA solution) 

0.47 < tan 2 $Q < 1.1 (LOW solution) 
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j 

1 (^ 1 ) 

!W 

!W 

i 

0.79 ±0.12 
-0.45 ±0.25 
0.34 ±0.29 

0.57 ±0.16 
0.49 ±0.28 

-0.60 ±0.23 

0.1 ±0.1 

0.69 ±0.18 

0.67 ±0.18 

1(e) 

2 (jj) 

3 (t) 


Table 1: Values of the MNS lepton flavour/neutrino mass mixing matrix U l3 as derived from solar 
and atmospheric neutrino oscillation data [2], 


Allowing also for the coupling of the state v 3 to muons, as suggested by data on oscillations 
of atmospheric neutrinos, but still setting U e3 = 0 [7] gives, in place of Eqn(8): 


(cos #i 2 - sin #i 2 ) cos 9 23 


(9) 


Assuming sin # 12 = 1 / 2 , and sin #23 = cos #23 = 1 /n/ 2 as suggested in Reference [7] 
(values consistent with the MNS matrix elements shown in Table 1) gives the value 2 
for the LHS of Eqn(9). It is clear, from these considerations, that the hypothesis that a 
coherent ‘lepton flavour eigenstate’ is produced in pion decay is experimentally excluded, 
with an enormous statistical significance, by the experimental measurements of R e /n and 
the MNS matrix elements. 

The correct prediction of the pion decay widths, that assumes incoherent production 
of the different neutrino mass eigenstates, is: 


T(vr —> l u) = T(vr —> £ u\) ± T(7t —> £ z/ 2 ) + r(7T —> £ u 3 ) 

* Ml(\Un\ 2 + \Un\ 2 + \U, 3 \ 2 ) = Ml (10) 


where the unitarity of the MNS matrix is invoked. The predicted value of R e / M , given by 
setting #12 = 0 in Eqn(7), is then independent of the values of the elements of the MNS 
matrix. Unlike in the case where the coherent production of a ‘lepton flavour eigenstate’ 
is assumed, no information about the MNS matrix can be obtained from measurements 
of .Rg/^ 

Giunti has claimed [8] that the argument just presented is flawed and the coherent 
states of massive neutrinos are produced in weak decay processes. To substantiate this 
claim it is suggested to define a ‘lepton flavour eigenstate’ not according to Eqn(3) but 
by instead writing the pion decay amplitude as: 


M.% = -M-iiUu ± + ■M.esUze. 
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fnm n V ud ^f,( 1 - 75 )^ 


( 11 ) 


where 


= UeiUitfo 1 + U^U2ti>u2 + U^ 3 U 3 ^ip V3 
— \Ue i\ 2 ' t Pv 1 + 11 T n | 2 ^v 2 + Wa\ 2 ^v 3 

^ (\Ua\ 2 + \Un\ 2 + \U(. 3 \ 2 )ipv 0 = (12) 

where r/’vo denotes the wavefunction of a massless anti-neutrino, the unitarity of the MNS 
matrix: \Un\ 2 + 1 1/^2 1 2 + \Ue 3 \ 2 = 1, has been used, and, as in Eqn(5), the purely kinemat- 
ical effects of non-vanishing anti-neutrino masses are neglected. Since the MNS elements 
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do not appear in Eqn(ll), the prediction given by this equation for R e /n is the same as 
the text-book massless anti-neutrino result, which is in excellent agreement with experi¬ 
ment and provides no information on the values of the MNS elements. However, since the 
amplitude (11) has no dependence on the values of the MNS elements, so that, unlike the 
correct Standard Model amplitude (2), the neutrino mass eigenstates are absent, it does 
not predict neutrino oscillations following pion decay. This seems now to be experimen¬ 
tally excluded by the observation of such oscillations in both atmospheric neutrinos [9] 
and the K2K [10] experiment. Actually, the anstaz of Eqn(ll) which seems to have been 
constructed precisely to avoid the constraint provided by Eqn(7), is in contradiction with 
the correct Standard Model expression for the pion decay amplitude, (2), which is lin¬ 
ear, not quadratic, in the MNS elements, and does contain the wavefunction of the mass 
eigenstate Xq. As shown below, neutrino oscillations are predicted by the amplitude (2), 
if eigenstates with non-degenerate masses exist. 

Although the incoherent nature of the production of the different neutrino mass eigen¬ 
states, as exemplified in Eqn(10) above, was pointed out more than 20 years ago by 
Shrock [11, 12], and the unphysical nature of coherent states of neutrinos of different 
mass was also discussed in the literature [13] the production of a coherent ‘lepton flavour 
eigenstate’ at a fixed time remains the basic assumption, in the literature, for the calcu¬ 
lation of the phase of neutrino oscillations [7]. The phase of each mass eigenstate is then 
assumed to evolve with time according to the Schrodinger Equation in its own rest frame 
(the so-called ‘plane-wave’ approximation). The assumption that all mass eigenstates are 
produced at the same time implicitly assumes equal velocities, since there is evidently a 
unique detection event at some well defined point in space-time. Still, in the derivation 
of the phase, the neutrino velocities, as defined by the kincmatical relation: v = p/E 
are assumed to be different. Thus contradictory hypotheses are made in space-time and 
in momentum space. These assumptions lead to the so-called ‘standard formula’ for the 
oscillation phase [7]: 

4>fr d = + 0(m 4 ) (13) 

J 2pu 

where rn l (i — 1, 2, 3) are the neutrino masses, L is the distance between the source pion 
and the detection event and is the neutrino momentum in pion decay to a massless 
neutrino. Units with h — c — 1 are used. 

It has been shown recently [14] that the factor of two in the denominator of this equa¬ 
tion is a consequence of a common production time for all mass eigenstates, that follows 
from the incorrect assumption that a coherent ‘lepton flavour eigenstate’ is produced. As 
will be demonstrated below, allowing for the possibility of different production times, clue 
to the incoherent nature of the neutrino production process, doubles the contribution of 
neutrino propagation to the oscillation phase. 

Other kinematical approximations that have been made in the literature, e.g. equal 
momenta and different energies or equal energies and different momenta, result in changes 
of only 0(m 4 ) to the phase, as compared to the exact result obtained by imposing both 
energy and momentum conservation in the neutrino production process. Thus all these 
results are equivalent at 0(m 2 ). 

In several previous papers written by the present author [14, 15, 16] results for the os¬ 
cillation phase differing significantly from Eqn(ll) have been obtained. These calculations 
were based on the covariant Feynman Path Integral formulation of quantum mechanics. 
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A general discussion of the foundations and applications of this approach is given in [17]. 
Here a concise re-derivation of the probability for e~ production resulting from the in¬ 
teraction of neutrinos from the decay n + —> /i + v at rest will be given. For simplicity, 
the effect of the finite pion lifetime as well as that of the smearing of the physical masses 
of the 7r and /j is neglected. As shown in Reference [15], the corresponding corrections 
to the oscillation formula are completely negligible. The basis of the calculation is very 
simple. Each interfering path amplitude has a common initial state (a pion that comes to 
rest in a stopping target at laboratory time t = to) and a common final state (defined by 
the neutrino detection event). There is one such path amplitude for each neutrino mass 
eigenstate. As neither the neutrino’s flavour nor its production time is observed, the cor¬ 
responding amplitudes must be coherently summed. Quantum mechanical coherence is a 
property of the path amplitudes (because they all have the same initial and final state) 
not of the neutrino production process. Indeed, the calculation is a direct application of 
the fundamental Eqn(7) of Feynman’s original Path Integral paper [18]. 

As in the usual ‘plane-wave’ derivation of Eqn(13) [7], the space-time propagators of 
the neutrino and the source pion can be simply obtained by solving the time-dependent 
Schrodinger Equation® in the particle rest frame. The path amplitudes for the detection 
of an e~ produced in the quasi elastic scattering process: un —» e~p by a neutrino from 
the decay at rest: tt + —> p + v are then: 

A efl (i) = A 0 U ei exp{-imi[T D (v) — r^z/)]}^ exp{—fm^fr^yr) — r 0 (7r)]} 

= AoUeiU^i exp(—iA0j) (* = 1,2,3) (14) 

where: 

A 0 =(e-p\T\un)(up + \T\7r + ) (15) 

The ‘reduced’ transition amplitudes (e~p\T\un ), (vp + \T\n + ) are defined by factorising 
out the MNS matrix elements, e.g.: 

(zy/i + |T|7r + ) = Up + \T\ii + ) (16) 

In Eqn(14), To(7r), Tiiy ) and Td{v) are the respective proper times of the n or v at 

which the pion comes to rest in the stopping target, the neutrino z/j is produced and the 
neutrino is detected. The overall propagator phase of the amplitude A efl (i) is the Lorentz 
invariant quantity: 


A 4>i = m\r D {v) - n(v)] + - r 0 ( n) 


(17) 


Since the pion is at rest: 


R.(vr) — T 0 (7r) — tj. — t 0 — tv — t 0 -— t D — t 0 

Vi 


1 + ^i- 
2 Pi. 


+ 0(m 4 


where ry is the velocity of z/j. The exact kinematical relation for the neutrinos: 

mL 


. At mL 
At = — =- 

7 Ev 
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(18) 


(19) 


“Note that this equation differs from that conventionally used in atomic physics by the inclusion of the rest-mass 
energy of the particle. Since the latter is at rest, the kinetic energy vanishes, and the equation is relativistically exact: 
sjm 2 + p 2 = m + T = to when p = T = 0. 


5 



and Eqn(18), substituted in Eqn(17), then give: 


A & 


mf 

Pu 


2 \p v 


L + m n (t D — t 0 — L) + 0(m 4 ) 


( 20 ) 


Eqns(14) and (20) then give, for the probability of e detection distance L from the 
source: 


Pen(L) — \A e n(l) + A eiJt (2) + A efJ ,(3)\ 

= 4|j4o| 2 U^U^sm^ + U^U^sin 

+u:,u:,u el u„ lS in 2 ^ 


2 023 
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L + 0(m 4 ) 


Am| = in 1 , - rrij 


( 21 ) 


( 22 ) 


(23) 


For neutrinos produced in the decay of an arbitary particle at rest [14, 15], or the 
decay, in flight, of an arbitary ultra-relativistic particle [15], Eqn(22) generalises to: 


'y 


AmU 1{- 

_y_ m 

Pu (1 - R 2 m) 


+ 0{m 4 


(24) 


where R m = rrix /% and rrix is the mass of the particle (or system of particles) recoiling 
against the neutrino produced in the decay of a source particle of mass mg. The neutrino 
is assumed to be light so that mg m*. For /3-decay at rest, the oscillation phase is 
given by [14, 15]: 


0y 


A ml 

Pu 


E/3 

Zpv 


L + 0(m 4 ) 


(25) 


where Ep is the total energy release of the /3-transition. 

The first term in the square brackets in Eqns(22) and (25) is the contribution to the 
oscillation phase of the neutrino propagator. It is factor of two larger than the similar 
contribution given by the standard formula (11). The second term in the square brackets 
in Eqns(22) and (25) is the contribution to the oscillation phase from the propagator 
of the source particle. It is opposite in sign to the neutrino contribution. It can be 
seen from Eqn(24), that, for light recoil masses such that R rn <C 1, the oscillation phase 
tends to vanish, strong cancellation occuring between the neutrino and source particle 
contributions. 

Note that the phase of the propagator of the source particle is the same, at the 
same time, in the different path amplitudes. The source may therefore be described as 
‘coherent’. Since, however, the decay process amplitudes are incoherent (see Eqn(10) 
above) the decay may occur at different times in different path amplitudes. It is this 
time difference that is at the origin of the contribution of the source particle propagator 
to the oscillation phase. In physical optics, the coherent source actually gives the entire 
interference phase, the contributions of photon propagators vanishing [17]. 
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Examination of Eqn(21) shows that the mechanism that governs the value of P e is 
interference between the path amplitudes for different neutrino flavours. A small value 
of P efl is not necessarily an indication of an approximate conservation of lepton flavour, 
but may be due to strong destructive interference between the different path amplitudes, 
independently of the values of the MNS matrix elements. This is made clear by consid¬ 
eration of the two flavour mixing scenario where U e 3 = 0 and Eqn(21) simplifies, using 
Eqn( 6 ), to: 


Pe^L) = |A 0 | 2 sin 2 26*12 sin 2 ^ = ^|A 0 | 2 sin 2 26* 12 (1 - cos(/q 2 ) (26) 

The term — cos 0i 2 in Eqn(26) originates in the interference of the path amplitudes corre¬ 
sponding to zq and z/ 2 . For small values of L, e~ production is suppressed by the almost 
complete destructive interference of these amplitudes, independently of the value of 0 V2 
i.e. of the degree of non-conservation of lepton number. The destructive nature of the 
interference is clue to the minus sign multiplying sin d 12 in the second row of the matrix on 
the RHS of Eqn( 6 ). This, in turn, is a consequence of the unitarity of the MNS matrix. 

Indeed, nowhere in the description of the ‘z/ e appearence’ experiment, just presented, 
do ‘lepton flavour eigenstates’ occur, although such an experiment is typically referred 
to [7] as ‘zq —» v e flavour oscillation’. In fact, only the mass eigenstates zq, zq and z/ 3 
appear in the amplitudes of the physical processes which interfere. It is the interfer¬ 
ence of these amplitudes in the production of the detection event that constitutes the 
phenomenon of ‘neutrino oscillations’. Indeed, no temporal oscillations of ‘lepton flavour’ 
actually occur. Still the terms ‘z/ e ’, ‘zq’ and ‘zq’ may still have a certain utility as mnemon¬ 
ics, even though they do not represent physical neutrino states for massive neutrinos. For 
example, it makes sense to refer to solar neutrinos, in a loose way, as a ‘z/ e beam ’ since 
the different physical components are all created together with an electron. Similarly, 
atmospheric neutrinos are predominantly ‘zq’, i.e., born together with a muon. 

The values of the correction factor C that must be applied to a value of Am 2 - obtained 
using the standard formula ( 11 ), to obtain the corresponding quantity as derived from the 
path amplitude formulae of Eqns(22), (24) or (25) are presented in Table 2 for different 
neutrino sources. The correction factor is: 

4>T d = (Am 2 ,)^ Es _ , = 1 - R 2 m 

4>ij Ami '^p v 2 74 

where Eg = Trig for particle decays and Eg = Eg for /3-decays. 

For muon decays and purely leptonic r decays, at the mean decay neutrino energy, 
as in the first row of Table 2, the correction to (A m ‘j.y tand is only 16 %. The small 
recoil masses in the decays: n —» eu, K —> eu lead to very large values of C (i.e. small 
values of the oscillation phase). For the electron capture process in the last row of 
Table 2, which is important for solar neutrino production, the correction factor is unity. 
The same is true in /3-decay processes when jq/p™ ax = 2/3, close to the mean value 
( p u ) = (11/16 )p™ ax of the spectrum of allowed /3-decays. There is then hope that mass 
differences derived from reactor neutrino experiments [19] using the formula (13) will not 
be too different from those given by the path amplitude formula (25). Actually, the phase 
should be evaluated separately for each /3-decay process contributing to the reactor flux, 
and a suitable weighted average performed. This is evidently a complex and difficult 
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Process 

Kinematical Condition 

C 

/X —> evv 

t —> evv 

T —> [1VV 

(Pu) = 7m e /20 

1.16 

7T —> jlV 

- 

0.372 

7T —> ev 

- 

37300 

K —* /iv 

- 

10.4 

K —> ev 

- 

466718 

K —> 7 TflV 

nix = m-K + m n ( C max ) 

1.52 

K —> nev 

nix = + rn e (C max ) 

5.71 

z N a -> z " 3 N a + e+ + V 

Pu/pT x = 2/3 
e.g. E u = 9.3 MeV 
for 8 B — Be* + e + + v 

1.0 

e~ + z N a N a + v 

e.g. e~ + 7 Be Li + v 

1.0 


Table 2: Values of the correction factor C relating the oscillation phase, as calculated using the 
standard formula, to that calculated using path amplitudes for various neutrino sources. 

undertaking. In contrast the corrections for atmospheric neutrinos (from 7r —> yuz/ and 
/x —» evv decays) and high energy solar neutrinos, essentially uniquely from 8 B /3-decay 
(as shown in the last-but-one row of Table 2), can be applied in a straightforward manner. 

Direct experimental discrimination between the standard formula for the oscillation 
phase (13) and the path amplitude formula (24) is possible in long-baseline terrestrial 
oscillation experiments by comparing ‘z/ M disappearence’ as observed in neutrino beams 
produced by either 7r —> fiv or K —> /j,v decays. Because of the relatively smaller recoil 
mass in the second process there is a strong relative suppression of the interference phase. 

For fixed neutrino momentum and source-detector distance, and considering only —> v T 

oscillations’: 

023(K) = 23 M = O.O358023(tt) (28) 

t'K 

The K2K experiment with a 250 km baseline has recently found indications of l v^ disap¬ 
pearence ’ for neutrinos of average energy 1.3 GeV from the process 7 t ^ /jv [10]. Using 
the standard formula (13) the event best fit to the data gives A m\ z = 2.8 x 10~ 3 eV 2 and 
9-23 — 7t/ 4 [10]. With L = 250km, Eqn(13) gives a phase 0 2 3 — 1-37 rad. Inserting these 
values into the standard two flavour L v P: disappearence’ formula [7]: 

P^{L) ~ 1 - sin 2 2d 23 sin 2 ^ (29) 

gives a suppression factor of 0.6 as compared to the no-oscillation (0 2 3 = 0) case. The 
standard formula (13) predicts exactly the same phase and suppression factor for neutri¬ 
nos of the same average momentum produced in K —> /j,v decays. However, if Eqns(24) 
and (28) are correct, 023 (K) = 49 mrad and a reduction of only 0.06 % of the ‘v p flux’ 
is predicted. This crucial test of the basic quantum mechanics of neutrino oscillations 
requires the installation of separated low energy charged kaon beams in the JHF-SK 
complex [ 20 ], or in similar long-linebase neutrino facilities. 




Another possible test is to search for an energy dependence of the oscillation phase 
for atmospheric neutrinos. Dne to the thickness of the Earth’s atomsphere, primary 
cosmic ray interactions that produce the pion (‘z/ M ’) and muon and l v e ') neutrino 
source particles occur at high altitude (~ 20km).For low neutrino energies, the flux of 
is roughly double that of ‘z/ e ’ since most muons decay before striking the surface 
of the Earth. However, for high neutrino energies (> lOGeV) most muons cross the 
atmosphere without decaying due to the relativistic time dilatation effect, so that only 
pion decays contribute appreciably to the neutrino flux which is thus predominantly (‘z^’. 
The oscillation phase for in-flight decays, 0l F , has been derived, using the Feynman path 
amplitude approach, in [15] with the result: 



R 2 Am-; 

TO _ P 

(1 - R 2 m ) cos 9 U E U 


J_\ Am?- L 

2 C ) cos 0 V E V 


(30) 


where 9 U is the angle between the directions of the neutrino and its source particle in 
the laboratory and E v is the laboratory neutrino energy. For n — » fiv decays the factor 
(1/2 C) has the value 1.34, whereas for /i —> evv decays the value of C corresponding 
to (p v ) quoted in Table 2 gives the value (1/2 C) = 0.431. Thus the oscillation phase 
is predicted by the path amplitude calculation to be larger at high values of E u where 
pion decays are the predominant neutrino source. This prediction may be tested using 
neutrino interaction data from the IceCube neutrino detector [21]. In a ten year exposure 
to atmospheric neutrinos more than 7 x 10 5 events are expected with muon energies E / 
greater than 100 GeV [22], For such events neutrino oscillations are strongly suppressed 
due to the ( p v )~ l factor in oscillation phase formulae. However, for lower energies, 10 
GeV < E^ < 100 GeV, where much larger numbers of events are expected, oscillations 
can be observed and the detected neutrinos are expected to originate almost entirely 
from pion decay. Thus the prediction of the path amplitude formula can be tested by 
comparing values of Am| 3 obtained at low energies from water Cerenkov detectors such 
as SuperKamiokande, where pion and muon decays contribute almost equally, with future 
IceCube results dominated by pion decay neutrinos. 

Acknowledgements 

I would like to thank D.Rapin for information concerning primary cosmic ray interac¬ 
tions and A.Rubbia for an interesting discussion of the possiblity to test the oscillation 
phase formula using atmospheric neutrino interactions. 


9 



References 

[1] Z.Maki, M.Nakagawa and S.Sakata, Prog. Theor. Phys. 28 870 (1962). 

[2] M.C.Gonzalez-Garcia and Y.Nir, Rev. Mod. Phys. 75 345 (2003). 

[3] N.Cabibbo, Phys. Rev. Lett. 10 531 (1963); 

M.Kobayashi and T.Maskawa, Prog. Theor. Phys. 49 652 (1973). 

[4] W.Marciano and A.Sirlin, Phys. Rev. Lett. 36 1425 (1976). 

[5] T.Goldman and W.Wilson, Phys. Rev. D15 709 (1977). 

[6] Review of Particle Properties , H.Hagiwara et al, Phys. Rev. D66 010001 (2002). 

[7] B.Kayser, Neutrino Physics as Explored by Flavour Change in Phys. Rev. D66 
010001 (2002) (Review of Particle Properties), and References therein. 

[8] C.Giunti, Coherence in Neutrino Oscillations, http://xxx.lanl.gov/abs/hep- 
ph/0302045. 

[9] T.Kajita and Y.Totsuka, Rev. Mod. Phys. 73 85 (2001). 

[10] M.H.Ahn et ai, Phys. Rev. Lett. 90 041801 (2003). 

[11] R.E.Shrock, Phys. Lett. B96 159 (1980). 

[12] R.E.Shrock, Phys. Rev. D24 1232 (1981); D24 1275 (1981). 

[13] C.Giunti, C.W.Kim and U.W.Lee, Phys. Rev. D45 2414 (1992). 

[14] J.H.Field, Eur. Phys. J. C30 305 (2003). 

[15] J.H.Field, Eur. Phys. J. C37 359 (2004). 

[16] J.H.Field, Quantum Interference Effects in the Detection Probability of Charged Lep¬ 
tons Produced in Charged Current Weak Interactions, http://xxx.lanl.gov/abs/hep- 
ph/0303152. 

[17] J.H.Field, Ann. Phys. (NY) 321 627 (2006). 

[18] R.P.Feynman, Rev. Mod. Phys. 20 367 (1948). 

[19] K.Eguchi et al, Phys. Rev. Lett. 90 021802 (2003). 

[20] Y.Itow et al, The JHF to Kamioka Project, KEK report 2001-4, 2001, 
http://xxx.lanl.gov/abs/hep-ex/0106019. 

[21] http://icecube.wisc.edu/ 

[22] M.C.Gonzalez-Garcia, F.Halzen and M.Maltoni, Phys.Rev. D71 093010 (2005). 


10 



